%6
x = [2,3,4];
y = [3,4,6];
k = (y(2)-y(1))/(x(1)-x(2));
syms t;
if k == Inf
    f1 = (y(1)+y(2))/2;
else
    k1 = 1/((y(2)-y(1))/(x(1)-x(2)));
    y1 = (y(1) + y(2))/2;
    x1 = (x(1) + x(2))/2;
    f1 = k1*(t - x1) + y1
end
k2 = 1/((y(3)-y(1))/(x(1)-x(3)));
y2 = (y(1) + y(3))/2;
x2 = (x(1) + x(3))/2;
f2= k2*(t - x2) + y2;
X = solve(f1==f2,t)
z1 = matlabFunction(f1);
Y = z1(X)
r = sqrt((x(1) - X)^2 +(y(1) - Y)^2);
X = double(X);
Y = double(Y);
ANG1 = angle((2-X)+i*(3-Y));
ANG2 = angle((4-X)+i*(6-Y));
m = ANG1:0.01:ANG2;
xx = X + r*cos(m);
yy= Y + r*sin(m);
plot(x(1),y(1),'o',x(2),y(2),'o',x(3),y(3),'o')
hold on
plot(xx,yy)
axis equal



